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CHAPTER  I:  nmtODDCTZOH 


In  this  study  we  investigate  the  coupling  of  electromagnetic  energy 
into  a  dielectric  slab  waveguide  on  one  interface  of  which  a  perfectly 
electrically  conducting  (p.e.c.)  half  plane  resides.  The  geometry  is 
shown  in  Figure  1.1.  The  top  face  of  the  slab  lies  in  the  y*  B  0  plane, 
with  the  p.e.c.  screen  residing  on  this  face  in  the  region  x'  >  0.  The 
slab  has  a  thickness  t'  and  constituitive  parameters  w*>ere 

e 2 '  “«y  be  complex.  It  is  imbedded  in  a  lossless  medium  with  parameters 
(uq'.Ei');  both  media  are  homogeneous  and  isotropic.  A  plane  wave 
polarized  transverse  magnetic  (TM)  to  the  z'  direction  impinges  on  the 
top  face  of  the  slab  at  an  angle  9  (0  <  0  <  x).  We  assume  there  is  no 
variation  in  the  z'  direction,  so  the  problem  is  strictly  two- 
dimensional. 

The  motivation  for  studying  this  problem  is  two-fold.  First,  this 
study  is  a  continuation  and  extension  of  the  work  of  Coblin  (1983). 
Second,  the  geometry  considered  provides  a  canonical  problem  which  can 
be  useful  in  modeling  metal-dielectric  junctions.  With  a  solution  to 
this  problem  one  is  able  to  consider  more  complicated  structures 
embodying,  for  example,  discontinuities  in  the  slab. 

This  investigation  builds  on  the  contributions  of  several  authors. 
Heitman  and  van  den  Berg  (1975)  have  studied  the  half  plane  in  multi- 
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layered  media,  while  Cobiin  (1983)  has  solved  the  problem  of  exterior 
scattering  from  a  single,  lossless  slab  screened  by  a  half  plane.  Tamir 
and  Felsen  (1965a)  have  developed  the  asymptotic  analysis  of  a  line 
source  in  a  slab.  In  addition,  many  authors  have  studied  the 
propagation  of  electromagnetic  fields  in  dielectric  waveguides  (e.g., 
Harrington  (1961),  Hessel  (1969a)). 

To  solve  this  problem,  one  formulates  an  integral  equation  for  the 
current  on  the  p.e.c.  screen  by  enforcing  the  condition  that  the  total 
tangential  electric  field  be  zero  on  the  p.e.c.  screen.  This  integral 
equation  is  amenable  to  the  Wiener-Hopf  technique,  which  provides  a 
solution  for  the  Fourier  transformation  of  the  current  on  the  screen. 
Multiplying  this  current  by  the  transformed  Green's  function  and  then 
taking  the  inverse  Fourier  transformation  gives  an  expression  for  the 
scattered  electric  field  in  the  form  of  a  spectral  integral.  The 
solution  was  carried  to  this  stage  by  Coblin  (1983),  and  his  analysis  is 
repeated  here  in  order  to  establish  consistent  notation  and  to  provide 
completeness . 

Coblin  carried  out  an  asymptotic  analysis  for  the  fields  in  the  two 
half  spaces  exterior  to  the  slab  but  not  for  the  fields  in  the  slab 
itself.  The  purpose  of  this  work  is  to  complete  the  asymptotic  analysis 
and  provide  an  asymptotic  representation  for  the  fields  within  the  slab. 

A  detailed  description  of  the  electric  field  within  the  slab  is 
developed,  and  an  effort  is  made  to  determine  how  various  factors 
influence  this  field.  To  this  end  the  representation  is  developed  in 


3 

terms  of  contributions  from  the  field  that  would  be  present  in  the 
absence  of  the  p.e.c.  screen,  a  collective  lateral  wave  (steepest 
descent)  field,  a  geometrical  optics  field,  and  a  guided  wave  field. 
The  guided  wave  field  is  represented  by  a  modal  expansion,  from  which  a 
launching  coefficient  formulation  is  developed. 


CHAPTER  II:  MATHEMATICAL  DEVELOPMENT 


2.1  Preliminary  Developments 


2.1.1  Definitions  and  Conventions 


In  this  chapter  we  develop  the  solution  to  the  electric  field  in 
the  slab.  Before  beginning  this  development  we  must  first  state  the 

e 

conventions  and  definitions  which  are  germane  to  the  solution  procedure; 
every  attempt  is  made  to  remain  consistent  with  the  notation  of  Coblin 
(1983).  The  anaylsis  is  conducted  for  time-harmonic  fields;  a  time 
convention  of  e^Ut  is  assumed  and  suppressed  throughout.  With  this 
convention  Maxwell's  equations  have  the  form 


■f  4  4  ,  .  -4 

V  x  e  —  -m  -  juy  h 

V  •  h  =  0 

*  4  -t  ♦ 

V  x  h  =  j  +  jue  e 


(2.1) 

(2.2) 

(2.3) 

(2.4) 


where  e  and  h  are  the  electric  and  magnetic  field  intensities, 
respectively,  P  is  the  electric  charge  density,  j  is  the  electric 
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current  density,  m  is  the  nsgnetic  current  density,  and  m'  and  t'  are 
the  constituitive  parameters  of  the  medium. 

The  media  are  passive,  homogeneous,  isotropic,  and  nonmagnetic,  but 
the  slab  may  be  lossy.  Each  medium  can  be  characterized  by  a  wavenumber 
k* ,  where 


k'  =  u/m  '  e '  -  k'  -  jk! 

r  1 

and 

k*  >  0  ,  k!  >  0  . 
r  l  — 

The  analysis  is  carried  out  in  the  Fourier  transformation,  or 
spectral,  domain,  and  throughout  the  chapter  terms  given  in  lowercase 
letters  denote  spatial  domain  quantities  while  uppercase  letters  denote 
spectral  domain  quantities.  We  define  the  Fourier  transformation  pair 

00 

A(k^)  =  F{ a(x)  }  ■  ^7  [  e^k*X  d*  (2.5) 

1«+  jc 

A(k  )  e'JV  dk  (2.6) 

-«.+  jc  X 


where  c  is  a  real  constant. 
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2.1.2  Normalization 


For  the  sake  of  convenience  we  introduce  a  normalization  in  the 
spectral  domain.  Throughout  the  context,  primed  quantities  (e.g.,  x') 
denote  unnormal ized  quantities,  and  normalized  quantities  are  unprimed 
(e.g.,  k2).  Defining  the  free  space  wavenumber  k^  =  u/ MqEq  =  2ir/AQ, 


k'.  =  k  */(!.*£  .  ,  ut  =  nAw., 
3  0  J  J  3  0  3 


Ej  =  eoej*  j  =  1>2> 


Thus  , 


where 


k'.  =  k  *k . 
J  0  j 


k.  =  /u  .  •  e  .  . 
J  J  J 


In  order  to  be  dimensionally  consistent  we  must  normalize  the 
spatial  variables  to  the  free  space  wavenumber  as  follows: 


(normalized  distance)  =  k^  *  (unnormalized  distance). 


x  =  kQ*x'  . 


etc.  For  example,  if  the  slab  had  a  thickness  of  one  half  the  free 
space  wavelength,  the  normalized  thickness  would  be 


/ 


t 


V*’ 


This  normalization  ensures  that  k*x  =  k'*x' ,  etc.,  and  that  the 
normalized  quantities  are  dimensionless. 


2.2  Development  of  Integral  Equation 

2.2.1  Integral  Equation  for  Current  on  Slab 

« 

Boundary  conditions  require  that  the  total  tangential  electric 
field  be  zero  along  the  p.e.c.  half  plane.  For  the  TM  case  this  means 
that  the  total  z'-directed  electric  field  component,  e  (x* ,y'),  be  zero 
for  y'=0,  x'>0.  We  partition  e^Cx'.y')  into  the  sum  of  two  terms, 


e  ( x'  ,y '  )  =  e^Cx'.y')  +  eR(x',y')  , 
z  z  z 


(2.7) 


where  the  following  definitions  have  been  made: 


NP/  ,  , i 

e  ( x  ,  y  ) 
z 


[  the  total  field  which  would  result 
J  with  no  p.e.c.  screen — including  the 
)  incident  field  and  the  field  scat- 
[  tered  by  the  slab 


,(x’,y’)  = 


the  field  radiated  by  the  current 
on  the  p.e.c.  screen 


After  applying  the  normalization  discussed  in  the  previous  section. 


2  1 


k  =  k.sina  , 
x  1 


a  =  u  +  jv , 


(2.47) 


dk  =  k.cosa  da. 
x  1 


This  a  plane  is  periodic  along  the  real  axis  with  period  2”,  so  we 
are  only  interested  in  the  interval  u  e  C — tt  ,  ti  ) .  The  mapping  has  the 
effect  of  "unfolding"  the  branch  cuts  associated  with  ±k^  so  that  the 
bottom  sheet  of  the  spectral  plane  is  exposed  in  the  one-sheeted  a 
plane.  The  integral  may  have  singularities  arising  from  GCk^.y) /G+ (kx) 
on  the  now-exposed  bottom  sheet  in  the  form  of  leaky  wave  poles  in 
addition  to  the  geometrical  optics  pole  associated  with  (k^  +  kjCos6) 
and  surface  wave  poles,  also  associated  with  G(k^ ,y) /G+ (kx> .  Under  the 
transformation  (2.47)  the  integral  in  (2.46)  may  be  written  as 


eXx.y)  =  5 


r  •  \  -jk.xsina 

I  cosa  G(kj8ina,y)e  1 

G+(kjSina) ( sina  +  cos9) 


da 


(2.48) 


where  we  have  defined  the  coefficient  £  as 


5  = 


-jEg8^(k^co8e)G(k^cos0,O) 

2wG_(-kjCOS0) 


(2.49) 


The  path  T  is  shown  in  Figure  2.3  for  the  case  of  a  lossless  slab  (k^ 
real)  to  emphasize  the  location  of  r  with  respect  to  singularities  in 


the  a  plane. 
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E"<V’,)  -  „0/I?-J,(kx)G(kx.y) 

-jt  B  (k^co80)G(k^cos0 ,0)G(kx,y) 
/5»G_(!k1co8e)G+(kx)  (k^  +  k^ 00807“ 


The  radiated  field  is  found  by  taking  the  inverse  Fourier  transformation 
of  this  expression,  or 


eV.y)  =  F_1{Ex(x,y)}  *  ^-\  E^(kx ,y ) e'JVdkx 

‘  -CB+  JC 

-jE  8  (k  co80)G(k. CO80 ,0)  ,  G(k  ,y)e  ^ jkxx 

=  — '  G_(-klCos0l -  !pG+^kxifkx”'k1lo8eidk*  (2‘46) 

where  the  path  P  is  the  collapsed  strip  of  analyticity  shown  in  Figure 

2.2. 


2.3  Asymptotic  Evaluation  of  the  Radiation  Integral 

2.3.1  Mapping  into  the  Angular  Spectral  Plane 

For  the  purpose  of  asymptotically  evaluating  the  integral  in 
equation  (2.46)  via  the  steepest  descent  method,  we  map  the  spectral 
plane  into  an  angular  spectral  plane,  the  a  plane,  using  the 


trans  formation 
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B  (k  ) 

P(kx)  =  VVW  -  MV =  S-(V  -  rorr  (2-43) 

for  all  k  . 

x 

Applying  the  edge  condition  at  the  edge  of  the  p.e.c.  screen  and 

then  investigating  the  asymptotic  behavior  of  the  functions  in  (2.43) 

shows  that  P(k  )-K)  as  Ik  |-*<».  Therefore,  by  Liouville's  Lemma  (Hille, 
xx 

1959a), 


P(k  )  =  0 
x 


for  all  k 


This  allows  us  to  solve  for  J+(k^): 


S+(V 
MV  =  G+(kx) 


-j«0  (k^cos0)G(k1cose,O) 
MQ2irG_(-k^cose)G+(kx)  (kx  +  k^cose) 


(2.44) 


2.2.5  Radiation  Integral 

Having  found  expressions  for  the  Fourier  transformations  of  the 

current  on  the  screen  and  the  Green's  function,  we  can  find  the 

R 

transformation  of  the  radiated  field,  E  (k  ,y) ,  which  is  given  by 

Z  X 
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A  (k  ) 

5TTT-W  **-<*,>.  !  (2-39) 

X 

The  decomposition  of  Coblin  (1983)  is  valid  for  complex  kj,  and  can  be 
accomplished  by  inspection.  The  functions  S+(kx)  and  S  (k^)  are  given 
by 

A+(k  )  A  (k  ) 

S-(kx)  =  G  (k  )  “  G  (-k  cose)’ 

-  x  1 

A  (k  ) 

W  G_(-k1cos6) ' 

With  the  decomposition  in  hand,  we  can  rewrite  the  Wiener-Hopf 
equation  as 

B_(k  ) 

•W'W  -  W  =  S-(V  -  gJTtt*  (2-42) 

Note  that  the  left  hand  side  is  composed  of  functions  analytic  for 

t  >  t_  and  the  right  hand  side  is  composed  of  functions  analytic  for 

t  <  t+.  Therefore,  we  can  analytically  continue  (2.42)  into  the  entire 

spectral  plane  where  each  side  of  (2.42)  is  equal  to  an  entire  function 

P(k  ): 
x 


(2.40) 

(2.41) 


1  7 


I  1  /  J  o 

f-2/Zi V  8 


l-2/2ir  "1  jH(k  )  „ 

[-V-  7T  -  lc'G(kx’0)*  x  • 
lx 


G+<kx>  = 


i/J.H(kx) 


(2.34) 


\  /  iun 

~u  e  y 

2/ZifJ  /k1  - 


Re(kx)<0 


G  (k  )  =  G  (-k  ) , 
-  x  +  x 


(2.35) 


where 


H(k  )  = 

X 


Ln{FU(On 

coa  — J  +  k  sin 
2]  x 


dc 


(2.36) 


and 


6(C) 


j 


(2.37) 


The  validity  of  the  factorization  has  been  checked  numerically  for 
tj'  =  — j 1 ) e q .  The  r.m.s.  error  in  (G  -  G+G_)/G  was  observed  to  be 

_  5 

less  than  10  over  both  the  top  and  bottom  sheets  of  the  kx  plane. 

This  factorization  allows  us  to  rewrite  the  Wiener-Hopf  equation  as 

A  (k  )  B  (k  ) 

•WW  *  c-nrr  -  rnrr-  ’-<'<%•  <2-38> 

-  x  x 


We  must  now  decompose  the  quotient  A+/G_  into  two  functions  S+(kx) 
and  S_(kx)  such  that  S+(kx)  is  analytic  for  t  >  t_,  S_(kx)  is  analytic 
for  t  <  ,  and 
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F(k  )  =  26 

X 


(B1  “  e2)e  jB 2*  "  (61  +  62)e+jB?t 
^  -  B^^V-  (Bj  +  62)2e+jV 


(2.32) 


The  first  tern  of  GCk^.O)  factors  by  inspection  into 


As  lkxl-~°.  V^7  for  k2*  real  or  co«plex.  Since, 

on  the  top  sheet  of  the  spectral  plane,  Im(B2)  <  0,  Re(jB2t)~»  as 
Ik  I  —  .  Therefore, 


F(kx)  >r  26 


0  -  (B1  ♦  B2)e+-’62t 

!0  -  (B1  +  B2)2e+iB2t 


*  1. 


These  arguments,  put  forth  by  Coblin  (1983).  are  restricted  to  lossless 
media  and  are  made  for  kx  on  the  top  sheet  only;  they  are  unaffected  if 
the  slab  material  is  lossy.  The  leaky  wave  poles  for  the  slab  lie  in 
quadrants  two  and  four  of  the  bottom  sheet,  just  as  the  proper  surface 
wave  poles  lie  in  quadrants  two  and  four  of  the  top  sheet.  Therefore, 
the  factorization  is  valid  on  the  bottom  sheet,  also,  and  it  is  given  by 


1  5 


b_(x)  «  CejklX. 

Con*' ',mring  an  infinitesimal  amount  of  loss  in  the  exterior  Medium, 
b  (x)  must  go  to  zero  as  x-»-“.  Using  this  knowledge  we  determine  that 
B  (k^)  is  analytic  for  t  <  Is^-kj).  (See  Figure  2.1.) 

Since  we  have  already  shown  that  GCk^.O)  is  analytic  and  nonzero  in 

the  strip  Im^^)  <  t  <  Im(-kj),  the  Wiener-Hopf  equation  (2.30)  exists 

in  the  strip  ImC-kjCose)  <  t  <  Im(-k^).  When  the  infinitesimal  lose  in 
the  exterior  nedium  goes  to  zero,  this  strip  collapses  to  the  real  axis 
as  shown  in  Figure  2.2. 

The  next  step  in  the  Wiener-Hopf  solution  is  to  factor  G(kx,0)  into 

two  functions  G  (k  )  and  G  (k  )  such  that  G  (k  )  is  analytic  and  nonzero 

+  X  X  +  X 

for  r  >  x_.  G_(kx)  is  analytic  and  nonzero  for  t  <  t+,  and 
G(k  .0)  *  G  (k  )G  (k  ). 

X  +  x  -  X 

Coblin  (1983)  has  found  this  factorization  for  a  lossless  slab.  In 
order  for  this  factorization  to  be  valid  for  a  lossy  slab,  G(kx,0)  must 
be  free  of  zeros  and  singularities  in  the  strip  of  analyticity, 
t_  <  t  <  x+  (which,  as  we  have  shown,  it  is),  and  we  must  be  able  to 
separate  it  into  a  product  of  two  terms,  one  of  which  must  limit  to  the 
value  1  as  Ik^l-*00  in  the  strip.  We  can  write  Gtk^.O)  as 

0(k,-0>  ■  <2-31) 


where 
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This  allows  us  Co  rewrite  (2.12)  as 

00 

a+(x)  -  b_(x)  *  |  j+(xQ) g( xQ ,x; 0)dxQ , 

*  -^OO 

-*»  <  x  <  <”.  (2.29) 

The  right  hand  side  of  equation  (2.29)  is  simply  the  convolution  of  j+ 
and  g(  so  after  applying  the  Fourier  transformation  to  both  sides  of  the 
equation  we  obtain  the  Wiener-Hopf  equation 

A  (k  )  -  B  (k  )  =  J  (k  )G(k  ,0).  t  <  r  <  r  .  (2.30) 

+  x  —  X  «f  X  X  —  + 

• 

For  this  transformation  to  be  valid  the  functions  A+,  B_,  J+,  and  G  must 
be  analytic  in  the  strip  t  <  t  <  t+,  and  G  must  be  nonzero  in  the  same 
strip.  We  have  already  shown  that  A+  is  analytic  for  t  >  In^-kjCosQ). 

The  current  j(x,0)  must  have  the  same  asymptotic  behavior  as  x-**0  as  the 
scattered  field,  and  thus  it  must  have  the  same  asymptotic  behavior  as 
the  incident  field  in  order  to  satisfy  boundary  conditions  everywhere 
along  the  screen.  Thus  j(x,0)  has  the  same  asymptotic  behavior  as  x-*® 
as  a+(x) ,  and  J+(kx)  must  be  analytic  in  the  same  region  in  which  A+(k^) 
is  analytic,  namely,  t  >  Im(-k^cose). 

We  next  note  that  b  (x)  is  proportional  to  the  scattered  field  on 
the  top  face  of  the  slab  (y=0)  for  x  <  0.  Along  the  slab,  as  x*-“>  we 
expect  this  field  to  show  the  same  behavior  as  a  negative- tr  avel  ing 


plane  wave,  or 


1  3 


(Coblin,  1983).  This  will  be  discussed  after  we  evaluate  the  field 
integral . 

Considering  a  lossy  slab  and  introducing  an  infinitesimal  amount  of 
loss  in  the  exterior  medium,  we  see  that  either 
IsK-k^)  <  Im(pole8,  zeros)  <  ImC-kj) 
or 

ImCk^)  <  ImCpoles,  zeros)  <  Im(kj). 

Thus,  GCk^.y)  is  analytic  and  nonzero  for 
ImCkj)  <  t  <  Im(-kj). 


2. 2. A  Wiener-Hopf  Solution  for  Jz(kx) 


In  the  Wiener-Hopf  procedure  (c.  f,  Mittra  and  Lee,  1971a),  we  seek 
the  Fourier  transformation  of  equation  (2.12).  In  order  to  accomplish 
this  transformation,  we  must  extend  the  limits  of  the  integral  in  (2.12) 
to  include  all  x  (i.e.,  x  t  (-<=,®)).  Therefore,  we  are  forced  to 
introduce  the  functions  b_(x)  and  j  +  (x)  defined  by 


( 

b_  (  x)  =  j 

t 

j  o, 

i  unknown 

func  tion , 

x  >  0 

x  <  0  , 

( 

j  +  (x)  = 

'  jz(x). 

,  o. 

x  >  0 

x  <  0  . 

(2.28) 


correspond  to  surface  waves  for  a  dielectric  slab,  and,  for  y  =  0,  the 
zeros  correspond  to  surface  waves  for  a  grounded  dielectric  slab 


such  that 


a(x,y)  =  z  a  (x,y) 


e^Cx.y)  =  -jua^x.y) 


h(x,y)  =  — ~  V  «  (za^Cx.y)]. 


The  potential  a^(x,y)  is  a  solution  to  the  equation 


(  „i  o'1 

|  V  +  k.  |  a  ( x ,y )  = 

l  *y  i  }  z 


f  -6(x) 6(y)  ,  i  =  2 


0, 


i  =  1 


subject  to  the  boundary  conditions  vdriich  require  t^iat  a^  and  i 
derivatives  be  continuous  at  y  =  0  and  y  =  -t.  Application  of 
Fourier  transformation  to  this  equation  yields 

r  3j  2  2  1  I  _7T7  6(y)*  1  =  2 

I  +  (k2  -  kZ)|  A  (k  ,y )  = 

(  6y’  1  X  ;  Z  X  '  I 

l  0,  i  =  1 

and  the  transformed  Green’s  function  is  given  by 


G(k  ,y)  =  -juA  (k  ,y). 

X  Z  X 


Using  the  approach  outlined  by  Felsen  and  Marcuvitz  (1973a) 


1  1 

(2.16) 

(2.17) 

(2.18) 

(2.19) 

ts 

the 

(2.20) 

(2.21) 

,  we 


obtain 


also  in  the  region  in  the  kx  plane  in  which  this  transformation  is 
analytic.  Ihe  transformation  of  a+(x)  is  given  by 


A+(kx)  =  F{a+(x) }  =  N"P(x.O)  eJV  dx 

j2E06i(kiCose)G(kic°se.O)  x 

uTuT^TTk  +  k.cosei 

0  x  1 

(  1  -  ej(kx  +  kicose)b  )  (2.14) 

The  limit  in  (2.14)  exists  and  is  0  only  when 

t  >  Im(-k1cos0)  =  t_,  kx  =  a  *  jt. 

Thus,  the  analytic  function  A+(kx)  is  given  by 

j2E06i(kxcose)G(kxcose,0) 

A+^kx^  n  u/S”  (k  +  k  cose)  ’ 

0  x  1 

r  >  Im(-kjCO80).  (2.15) 


2.2.3  Fourier  Transformation  of  the  Green's  Function 

In  order  to  solve  for  the  current  on  the  half  plane,  we  must  find 
the  Fourier  transformation  of  the  Green's  function  for  a  line  source  on 
the  top  face  of  the  slab,  G(kx>y),  and  find  the  region  in  the  spectral 
plane  where  G(kx,y)  is  analytic  and  nonzero.  To  this  end,  we 
temporarily  introduce  the  concept  of  a  magnetic  vector  potential 
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ez(x' ,y')  =  uQ  j  jz(x0)g(x* ,x0;y' )  dXg. 


(2.11) 


Application  of  the  boundary  conditions  produces  the  integral 
equation 


NP  t 

ez  (x',0)  =  -uQ  I  jz(x0)g(x' ,xQ;0)  dxQ,  (2.12) 

0  <  x'  <  =o. 

This  equation  is  of  a  form  suitable  for  application  of  the  Wiener-Hopf 

« 

technique,  which,  when  applied,  will  give  a  solution  for  the  Fourier 
transformation  of  the  current  on  the  half  plane.  To  find  this  solution 
we  must  find  the  Fourier  transformation  of  the  no-plane  field  and  of  the 
Green's  function. 


NP 

2.2.2  Fourier  transformation  of  e  (x,y) 


To  apply  the  Wiener-Hopf  technique,  we  define  the  function 


a+(x) 


— e  (x,0)  ,  x  >  0 
4q/2tt  z 

0  ,  x  <  0. 


(2.13) 


We  are  not  only  interested  in  the  Fourier  transformation  of  a+(x),  but 


8 


NP 

the  e ,  or  "no-plane,"  component  of  the  field  can  be  found  by  solving 
the  problem  of  a  plane  wave  of  magnitude  Eg,  given  by 

inc,  v  _  _  jk,(xcoe0  ♦  yeinQ)  . 

e  (x,y)  =  E  eJ  1  3  ,  (2.8; 


incident  upon  a  dielectric  slab.  The  no-plane  component,  given  by 
Coblin  (1983),  is 


NP. 

e_  (x,0) 


E0B1(k1cose)G(k1co8et0)  ejklXC°8e,  (2.9) 


where  we  have  rearranged  the  expression  to  include  G(k  ,y),  the  Fourier 

•  x 

transformation  of  the  Green's  function  for  a  line  source  on  the  top  face 
of  the  slab,  and 

B . (k  )  =  / k?  -  k2  ,  i  =  1,2,  (2.10) 

1  x  IX 

where  we  choose 

Im(B^ )  <  0,  or,  if  Im(6  x )  =  0,  Re(B  ^ )  >  0. 


Choosing  the  sign  of  the  square  root  in  this  manner  ensures  that  the 
fields  satisfy  a  radiation  condition  at  infinity  and  defines  the  top 
sheet  of  the  spectral  plane. 

We  can  determine  e  from  the  convolution  of  the  current  induced  on 
z 

the  half  plane  and  the  Green's  function  for  a  z'-directed  line  source  on 
the  slab,  i.e.. 
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2.3.2  Deformation  of  V  to  the  Steepest  Descent  Path 

i 

To  evaluate  (2.48)  using  the  steepest  descent  method,  we  must  find 
the  path  along  which  the  exponential  in  the  integrand  of  (2.48)  decays 
most  rapidly — the  steepest  descent  path,  or  SDP,  which  passes  through 
the  saddle  point,  a^,  where  the  derivative  of  the  argument  of  the 
exponential  with  respect  to  a  is  zero;  i.e., 

•^(- jkjXsina)  =  0. 

a 

s 

The  saddle  points  in  the  region  of  interest  are  ag  =  n/2  for  x  >  0  and 
ag  =  -tr/2  for  x  <  0.  The  SDP  lies  on  the  contour  through  the  saddle 
point  characterized  by  the  imaginary  part  of  the  argument  of  the 
exponential  in  (2.48)  being  constant.  It  is  given  by  points  that 
satisfy 

Im(- jkjxsina)  =  Im(- jk^xsina^ )  =  -kj|x|  ,  (2.50) 

or 

sinh(v)  =  cot(u)  (2.51) 

where 

’  (0,n)  ,  x  >  0,  a  =  n/2 
s 

u  e 

(-n,  0),  x  <  0,  ag  =  -tt/2. 

Thus,  we  will  be  concerned  with  two  steepest  descent  paths,  SDP  through 
ag  =  -rr/2  and  SDP  through  ag  =  -w/2.  These  paths  are  shown  in  Figure 


2.4.  When  the  path  r  is  deformed  to  the  SDP,  singularities  will  be 

I  swept,  thus  contributing  residue  terms  to  e  (x,y)  in  addition  to  the 

contribution  from  the  evaluation  of  the  integral  along  the  SDP.  This 
r 

shows  that  e^  (x,y)  is  made  up  of  components  from  four  contributors: 

1.  The  result  from  the  evaluation  of  (2.48)  along  the  SDP. 

2.  The  residue  obtained  from  the  capture  of  the  geometrical  optics 
pole. 

3.  The  residue  obtained  from  the  capture  of  any  number  of  surface 
wave  poles . 

4.  The  residue  obtained  from  the  capture  of  any  number  of  leaky 
wave  poles . 

Tamir  and  Felsen  (1965b)  have  found  that  the  contribution  to 

R  SDP 

ez  (x.y)  from  the  SDP  integration,  (x,y),  is  equivalent  to  a 

collective  lateral  wave  as  shown  in  Figure  2.5.  The  contribution  from 

GO 

the  geometrical  optics  pole,  e^  (x.y),  accounts  for  the  shadowing  of 

the  incident  field  by  the  p.e.c.  screen.  The  contributions  from  the 

M 

surface  and/or  leaky  wave  poles,  e ^  (x.y),  can  be  characterized  as  a 
series  of  modes  propagating  and/or  decaying  in  the  slab  in  a  manner 
similar  to  that  in  a  dielectric  waveguide.  The  total  radiated  field  is 
found  by  summing  all  the  individual  contributors: 
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R,  N  SDP,  v  .  GO,  x  .  _  M,  x 

e  (x,y)  =  e  (x,y)  +  e  (x,y)  +  t  e  (x,y). 
z  z  z  z  ' 


Each  component  of  (x,y)  is  considered  in  the  following. 


2.3.3  SDP  Integral 


SDP 

The  contribution  from  e^  (x,y)  is  given  by  the  integral 


SDP,  -  r 

ez  (x.y)  =  ? 


,  .  x  -jk.xsina 

I  cosa  GCkj  sina,y)e  I 


J  SDP 


G+(k1sina) (sina  +  cose)  da’ 


This  integral  can  be  solved  in  a  straight-forward  manner  by 
steepest  descent  method  (Felsen  and  Marcuvitz,  1973b,  Mittra  and 
1971b).  We  note  that  (2.53)  is  of  the  form 


SDP,  \  x  xq( a) 

e  (x,y)  =  |  K(a,y)e  da 

Z  1  SDP 


where 


and 


cosa  G(k^sina.y) 

K(a,y)  r,  (k  sina)(sina  ♦  cose) 
+  1 


q(a)  =  tjkjSina  for  x  ^  0. 


(2.52) 


(2.53) 

the 
Lee , 

(2.54) 

(2.55) 


(2.56) 


Making  the  substitution 


s  =  q(cO  ~  q(<0 


(2.57) 


we  note  that  Im{q(a)}  is  constant  along  the  SDP,  so  s  is  real  along  the 
SDP.  This  allows  us  to  write  (2.53)  as 


SDP ,  ,  _ 

e  (x.y)  -  e 
z 


-jlkjxl  f 


|  M(s,y)e  XS  ds  , 


(2.58) 


M(s.y)  =  K(a,y)  ^ 


(2.59) 


a  on  SDP 


We  can  express  this  integral  as  an  asymptotic  series  for  large 


|kjx| : 


“jlk.xl  7  |M2n(0.y)[  -d"  f  t  |j 


SDP,  \  -  j  |  k,  x  |  _  M  (0,y)  -d  ir  V 

e  (x.y)  *  e  J  1  Z  TTYT - Z — I  ( • 

Z  n=0 j  2n [(dk  x)niklX;  Jf 

v  1  ) 


as  |k^x|-«=. 


(2.60) 


The  first  term  of  this  series  is  zero  due  to  the  factor  cos(a),  which  is 

zero  for  a  =  a  ,  so  the  second  term  is  the  dominant  one,  and 
s 


SDP  j  f’kl  -»/» 


.in  . 

+  J  {- j  -  lkix ! 


(2.61  ) 


where 
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Y(y) 


_  ; 


-H(k. ) 
-e  1 


2/ky(cos0  +  l)co8  C &2i 


•sin(e21y) , 


j4^eH(V 

1  621 ^CO60  "  1 ^ 


tanle  (|)]  t 

- r — *sin[621(yh5-)]  + 


cos[  B21  (-|)  ] 


cot[ 6,.  0^)] 

— — — •co8[a21(y4I)] 


sin[821(‘2>] 


x  >  0 


(2.62) 


x  <  0. 


and 


B21  "  62(kl)* 


(2.63) 


This  expression  represents  a  collective  lateral  wave  made  up  of  a 
finite  summation  of  lateral  waves  as  shown  in  Figure  2.5  (Tamir  and 
Felsen  1965).  Each  lateral  wave  in  the  collection  propagates  for  a 
unique  distance  in  the  exterior  medium  and  undergoes  a  unique  number  of 
reflections  into  the  dielectric  slab  at  the  critical  angle.  This 
critical  angle  is  simply  the  angle  of  total  internal  reflection  0^  given 


sin(e  )  =  /e./t,  . 
c  12 


by 
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2.3.4  Uniform  SDP  Result 


The  contribution  from  e^  (x,y)  exhibits  singular  behavior  for 
incidence  angles  near  *  for  x  >  0  and  for  incidence  angles  near  0  for 
x  <  0.  To  correct  for  this  behavior,  the  uniform  SDP  evaluation  method 
(Felsen  and  Marcuvitz,  1973c)  is  applied  and  leads  to  a  result  that  has 
no  singular  behavior  at  these  angles.  It  is  given  by 


[  SDP ,  '  ^88n(x)E0 

e  (x,y)  s  - 

{  z  J  Junrf  u 


■6^(k^cose)G(kjCose.y)x 


/r, - ; - : - — —  jn/4l  jk.xcose  • 

2Q | / 1 k^ x |  +  k^xcosS  e  eJ  1  + 

i. 

e~j|klx|  (  - jx/4  _  ej*M  \ 

/jk^xl  +  kjcose  [  6  2( | kjX |  *  kjXcose) j 


.  { 3u 


-j|—  +  |kjX  I  j 


(2.64) 


where  Y(y)and  are  given  in  (2.62)  and  (2.63),  respectively,  and 


Q(y)  =  j  e  dx. 

y 


(2.65) 
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2.3.5  Geometrical  Optics  Contribution 


Hie  geometrical  optics  field  accounts  for  the  shadowing  effect  of 

the  screen.  It  arises  from  a  singularity  in  the  integrand  of  (2.46), 

the  geometrical  optics  pole,  at  =  -k^cosG.  In  an  asymptotic  solution 

we  would  expect  to  encounter  this  field  when  x  >  0  and  not  vftien  x  <  0. 

This  is  exactly  what  we  find;  as  a  consequence  of  the  Wiener-Hopf 

solution  procedure,  the  geometrical  optics  pole,  a  =  6  -  ir/2,  always 

8  ® 

lies  below  the  path  T  and  in  the  interval  [-ir/2,it/2]  on  the  real  axis  of 


the  a  plane.  Thus,  it  is  swept  by  the  deformation  to  SDP  (i.e., 
observaton  points  under  the  half  plane)  and  not  by  the  deformation  to 
SDP  .  This  pole  is  encircled  in  the  negative  (clockwise)  direction, 
therefore  the  geometrical  optics  contribution  is  given  by 


e^°(x,y)  =  -j2u£U(x)*Re6 


_  .  \  -ik.xsina 

CosaG(kj8ina,y)e  J  1 

^  G+(kjSina) (sina  +  cose)  '  ago 


.  2/2T  _  „  .  „  v  jk.xcosG 

U(x)* -  E-B, (k.cosG)G(k  cosG,y)eJ  1 

u  0  1  1  1  J 


where  U(x)  is  the  unit 


U  ( x) 


1. 

0. 


step  function, 
x  >  0 
x  <  0. 


(2.66) 


(2.67) 


We  note  that 
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GO  , 
e,  (x.y) 


NP, 

-e,  (x.y) 


for  x  >  0.  This  is  as  expected,  and  shows  that  the  geometrical  optics 
field  does  indeed  account  for  the  shadowing  from  the  screen  by  canceling 
the  no-plane  field. 


2.3.6  Modal  Contributions 

« 

The  remaining  contributions  to  the  radiated  field  arise  from  poles 
of  G/G+  in  the  integrand  of  (2.46).  Both  surface  wave  poles  on  the  top 

sheet  and  leaky  wave  poles  on  the  bottom  sheet  may  be  swept  in  the 

.  +  -  . 

deformation  from  r  to  SDP  or  SDP  ,  each  pole  contributing  a  residue 

£ 

term  to  e^  (x.y).  Each  of  these  residue  terms  can  be  expressed  as  a 
product  of  a  launching  coefficient — L(©,Op),  a  transverse  (y  direction) 
modal  function — *(y,cO,  and  a  longitudinal  (x  direction)  propagation 
factor — an  exponential  function  of  x  and  a^,  where  denotes  the  value 
of  the  pole  in  the  a  plane.  With  this  representation  we  have  an 
efficient  form  to  study  separately  the  modal  structure  of  the  field,  the 
factors  that  influence  the  launching  of  a  mode,  and  the  propagation  of  a 
mode . 

For  each  pole  of  G/G+  swept  in  the  deformation  to  the  steepest 


descent  path. 
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M,  ^ 

ez(x,y) 


i  -  \  “jk.xsina 

cos  aG(k^ sina,y)e  1 


j2n5  Res  I  G+(kj8ina) (sina  +  cos0) 


(2.68) 


for  x  0. 


The  evaluation  of  the  residue  is  tedious  but  straightforward  and  the 
resulting  expression  can  be  rearranged  into  the  desired  representation. 


e”(x.y)  =  EqL( 0,a  )®(y,a  )e  ^kxX 


(2.69) 


where  kx  =  kjsina^.  The  modes  that  result  can  be  segregated  into 
readily  interpretable  classes.  Under  the  half  plane,  i.e.,  x'  >  0,  the 
electric  field  must  vanish  at  y'  =  0,  and  all  modes  are  governed  by  thi6 
condition.  For  x'  <0  the  slab  is  an  open  structure,  and  a  mode  set 
manifesting  both  odd  and  even  symmetry  about  the  center  of  the  slab 
(y*  =  -t'/2)  arise.  Specifically,  for  x'  >  0 


1.(0, a  )  = 


2  ~H ( ) 

8.  (k  cos  9)  cot(  a  )B„  G(k,cos0,O)e  x  ,  ~ 

11  p  2p  1  i'-2/2ii 


G  (-k. cos  0) (k, cos  0  +  k^)sin(6„  OvtT +  k^  1 
1  1  x  2p  1  x 

{lV°8(B2pt)  '  8lpB2PtSin(62pt)]  + 

j[Bjptcos(B2pt)  ♦  B^sinlB^t)]}'1, 


l  <*>  ) 


(2.70) 
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*(y,a  )  =  sin( 8,  y) , 

P  ZP 

where 

8.  *  8.  (k. sina  ) , 

lp  11  p 

8.  =  80(k  sina  ) . 

2p  2  1  p 

For  x'  <  0  odd  nodes  are  given  by 


(2.71) 


(2.72) 

(2.73) 


8.  (k.cos 6)cot(a  )B0  >4c,  -  k^G(k,  cos  6, 0)e  x  ,  r,m~ \/ 1 

. \  _  1  1  p  2p  1  x  1  J-2 /2 ir ! 

1.(9, a i  )  =  - - - c K - - - -  -  xA, 

P  G  (-k1co8  0)(k,cos0  +  kp)8in[B0  (^)]  u 

-  1  1  x  2p  / 

(2.74) 

A=  {2l(6lp  +  B2p)co8(62pt)  "  BlpB2ptsin(fV)]  + 

^4BlpB2pSin(  B2pt^  +  Blpt^Blp  +  B2p^c°6^  B2pt^  ^  ’  (2*75) 

#(y,a^)  =  sin[ 82p(y+j)J  (2.76) 


with 


lp 


and  8„  given  above, 
ZP 


while  the  even  modes  are  given  by 


L(9,a  ) 
P 


8. (k. cos  9) cot( a  ) 6„  /k, -  kPG(k, cos  9,0)e 
1  1 _ p  2p  1  x  1 _ _ 

G_(-k ^cos  0) (k^cos  9  +  kP)cos[ B2p(y)] 


-H(kP)  1/2 

X  A 

- X  ( - ^ - J  *A' 


*(y,a  ) 
p 


cosi  B2p(y+^H 


(2.77) 

(2.78) 


wi th  8 .  , 

lp 


80  ,  and 
2P 


A  given  above. 
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The  reader  is  reminded  that  these  are  large-x  asymptotic  field 
representations.  The  continuity  of  e^  at  x'  =  0  is  embodied  in  the 
integral  representation  (2.48),  which  constitutes  a  complete  expression 
for  the  field.  The  changing  of  asymptotic  representations  between 
x'  >  0  and  x'  <  0  does  not  contradict  continuity. 

The  remainder  of  this  study  will  be  devoted  to  an  investigation  of 
the  poles  of  G/G+  and  the  modal  fields  which  arise  from  them. 


j 

CHAPTER  III:  ANALYSIS  OF  RESULTS 

3.1  Preliminary  Developments 

In  a  uniform  dielectric  slab  waveguide  of  thickness  t  and 
permittivity  e^,  electromagnetic  energy  can  propagate  as  guided  waves  in 
the  form  of  modes.  These  modes  will  propagate  in  the  slab  in  the 
longitudinal  direction  and  will  radiate  energy  into  the  exterior  medium 
in  vdtich  the  slab  is  imbedded.  They  can  be  separated  into  two  types: 
modes  that  are  odd  with  respect  to  the  center  of  the  slab  and  modes  that 
are  even  with  respect  to  the  center  of  the  slab.  If  the  slab  is 
grounded  by  placing  a  p.e.c.  plane  on  one  surface,  image  theory  tells  us 
that  this  must  be  equivalent  to  a  dielectric  slab  waveguide  of  thickness 
2t,  but  since  each  mode  must  individually  satisfy  boundary  conditions 
along  the  ground  plane,  only  those  modes  that  are  even  with  respect  to 
the  ground  plane  may  exist  (Figure  3.1).  Both  types  of  behavior  should 
be  evident  in  the  problem  under  consideration.  For  large  |k^x  |  the 
portion  of  the  slab  under  the  p.e.c.  half  plane  (x  >  0)  should  behave  as 
a  grounded  dielectric  slab  waveguide,  and  the  portion  of  the  slab  not 
under  the  half  plane  (x  <  0)  should  behave  as  a  slab  waveguide  in 
uniform  space.  This  behavior  is  apparent  in  the  analysis  of  the  poles 
that  produce  the  guided  wave  terms  and  in  the  analysis  of  the  modal 
expressions  that  result  from  these  poles. 
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In  addition  to  the  geometrical  optics  pole,  two  types  of  poles  are 

encountered  in  the  evaluktion  of  e  ^(x,y).  Surface  wave  poles  (SWP's) 

z 

which  lie  on  the  top  sheet  of  the  spectral  plane  and  leaky  wave  pole6 
(LWP's)  which  lie  on  the  bottom  sheet  may  contribute  field  terms.  The 
two  types  of  poles  give  rise  to  two  different  types  of  waves  (Hessel, 
1969c).  A  surface  wave,  produced  by  a  SWP,  behaves  as  a  propagating 
mode  in  a  slab  waveguide.  It  propagates  in  the  longitudinal  direction, 
attenuated  only  by  losses  in  the  slab.  As  it  propagates,  it  radiates 
energy  in  the  exterior  medium.  This  radiation  decays  exponentially  away 
from  the  slab,  while  the  transverse  variation  is  sinusoidal  inside  the 
slab  (Figures  3.1  and  3. 2. a).  A  LWP  produces  a  leaky  wave,  which 
behaves  much  like  an  evanescent  mode  inside  the  slab  in  that  it  decays 
exponentially  in  the  longitudinal  direction.  Therefore,  it  does  not 
make  a  significant  contribution  to  the  total  field  in  an  asymptotic 
solution.  However,  since  a  LWP  occurs  on  the  bottom  sheet  of  the 
spectral  plane,  a  leaky  wave  will  not  satisfy  a  radiation  condition  in 
the  transverse  direction;  it  grows  exponentially  away  from  the  slab 
(Figure  3.2.b).  The  leaky  wave  field  will  satisfy  a  radiation 
condition,  though,  because  it  only  exists  in  a  restricted  portion  of  the 
spatial  plane  (Figure  3.2.c).  In  Figure  3.2.c  the  center  line  of  the 
6lab  and  the  direction  of  growth  of  the  leaky  wave  define  an  angle 
Og  ^  but  the  leaky  wave  pole  only  contributes  a  field  term  inside 

the  wedge  defined  by  the  center  line  and  the  angle  scontr  •  The  angle 
^growth  alway8  8reater  than  econtr  •  ensuring  that  the  leaky  wave 
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field  does  satisfy  a  radiation  condition  in  the  region  in  vrtiich  it 
exists. 
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map  into  the  shaded  region. 
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P 


for  the  radiation  integral  showing 
rip  of  analvticity. 


k|COs# 
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Hie  physical  properties  of  the  slab  and  the  location  of  the  guided  wave 
poles  also  had  an  effect  on  the  launching  of  the  modes.  The  large 
magnitudes  of  the  launching  coefficients  was  unexpected,  and  at  present 
an  adequate  explanation  is  not  available.  Further  study  of  this 
phenomenon  is  warranted. 

With  this  solution  completed  for  the  geometry  considered  in  this 
study,  one  is  prepared  to  extend  the  work  to  more  complicated  slab 
structures  or  to  delve  deeper  into  the  physical  mechanisms  at  work  in 
the  scattering  process. 


CHAPTER  IV:  CONCLUSION 

In  this  study  we  have  investigated  the  scattering  of  an 
electromagnetic  plane  wave  by  a  dielectric  slab  shielded  by  a  p.e.c. 
half  plane  residing  on  the  top  face  of  the  slab.  We  developed  solutions 
for  the  Green's  function  for  a  line  source  on  the  slab  and  for  the 
current  on  the  half  plane,  extending  the  results  to  accomodate  losses  in 

the  slab.  The  total  field  was  expressed  as  the  sum  o_f  a  no-plane  field, 

which  accounted  for  the  effects  of  the  slab  and  the  incident  plane  wave, 
and  a  radiated  field  which  accounted  for  the  effect  of  the  screen.  This 

radiated  field  was  developed  asymptotically  as  the  sum  of  a  collective 

lateral  wave  contribution,  a  geometrical  optics  contribution,  and 
contributions  from  guided  wave  poles.  The  contributions  from  the  guided 
wave  poles  were  expressed  as  modes  giving  rise  to  surface  waves  or  leaky 
waves .  We  were  able  to  state  these  modes  as  the  product  of  a  transverse 
modal  function,  a  propagation  factor,  and  a  launching  coefficient.  This 
launching  coefficient  played  a  similar  role  in  the  launching  of  these 
guided  waves  to  that  of  a  diffraction  coefficient  in  diffraction  theory. 
The  launching  of  the  guided  modes  was  investigated  in  some  detail,  and 
an  effort  was  made  to  explain  the  physical  mechanisms  that  affected  the 
launching  of  these  guided  modes.  It  was  observed  that  the  edge 
diffraction  effect  played  a  major  role  in  the  mode- launching  phenomenon. 
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Another  interesting  effect  is  seen  in  the  launching  coefficients 
for  a  slab  in  free  space  with  permittivity  10  £q  and  thickness  Aq/2 
(Figure  3.8).  For  this  case  the  launching  coefficient  is  not  zero  for 
6=0  and  6  =  it.  Normally  8^(k^cos6),  which  is  zero  for  6=0  and 
6  =  ir,  forces  the  entire  launching  coefficient  to  zero  at  these  angles. 
For  these  parameters,  however,  there  is  a  pole  of  G(±kj,0)  that  cancels 
the  zero  of  6^(±k^).  Thi6  phenomenon  will  occur  whenever 
t/e2  -  =  nir,  n  an  integer. 

When  this  relationship  is  satisfied,  the  electrical  thickness  of  the 
slab  is  such  that  the  no-plane  field  is  no  longer  zero  for  6  =  0  or 
6  =  n,  and  modes  can  be  excited. 


shapes  are  relatively  independent  of  the  location  of  the  poles  that 
produce  them.  This  further  supports  the  hypothesis  that  edge 
diffraction  dominates  the  launching  of  modes  and  also  shows  that  the 


properties  of  the  slab  influence  the  shape  of  L(0,  cO  more  than  the 
location  of  the  poles. 

The  launching  coefficients  have  only  one  or  two  local  maxima  over 
the  range  of  6.  The  number  of  maxima  remains  the  same  for  each  pole 
found  for  a  given  set  of  parameters.  For  slabs  that  are  electrically 
thin  (e.g.,  '  =  4eq,  t'  =  0.2X0>  Ej  '  =  eq),  L(e,a^)  has  a  single 

maximum,  and  as  the  electrical  thickness  of  the  slab  increases  the 
second  local  maximum  first  appears  as  a  bend  in  the  coefficient  and  then 
as  an  extremum.  One  can  develop  a  relation  for  the  location  of  these 
maxima  by  setting  the  6  derivative  of  L(0,(O  equal  to  zero,  but  the 

resulting  equation  is  complicated  and  yields  little  if  any  physical 
insight . 

Even  though  the  location  of  does  not  significantly  affect  the 

shape  of  1.(0, ct  ),  a  has  a  large  effect  on  its  magnitude.  Surprisingly, 
P  P 

the  magnitude  of  L(0,a  )  is  much  larger  for  the  higher  order  modes  in  a 

P 

given  slab  than  for  the  lower  order  modes.  For  a  waveguide  one  sees  the 
opposite  effect — the  lowest  order  modes  in  a  waveguide  carry  the  most 
energy.  Also,  one  would  expect  L(9,a  )  to  have  a  magnitude  of  order 
one,  when  in,  fact,  its  magnitude  can  be  much  larger,  especially  for  the 
higher  order  modes.  At  present  an  adequate  explanation  for  these 
effects  i6  not  available. 
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expression  similar  to  the  modal  launching  coefficient  and  compare  the 

SDP 

two.  This  expression,  which  we  call  L  (6),  has  three  forms.  These 
forms,  which  are  comparable  to  the  three  forms  of  L(e,a  ),are 


Figure  3.7  shows  L(6,a^)  and  L  ( 9)  as  functions  of  6  for  a 

representative  case.  We  note  that  except  near  the  geometrical  optics 
SDP 

pole,  where  L  ( 9)  is  singular,  the  two  functions  have  the  same  general 
shape.  They  differ  greatly  in  magnitude,  but  this  would  be  expected. 


These  results  support  our  claim  that  edge  diffraction  dominates  the 


mode- launching  effect. 

The  launching  coefficients  are  smooth,  slowly  varying  functions  of 
9.  For  a  particular  choice  of  slab  thickness  and  permittivity,  their 
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The  transverse  modal  functions  also  give  the  expected  results.  For 
both  surface  waves  and  leaky  waves  they  show  the  odd  or  even  behavior 
described  previously.  If  the  slab  is  loss,  ss  the  transverse  variation 
of  a  surface  wave  is  strictly  sinusoidal  and  remains  essentially 
sinusoidal  for  a  lossy  slab  (Figures  3.3,  3.4,  and  3.5).  Leaky  waves 
exhibit  an  overall  growth  in  magnitude  in  the  transverse  direction  away 
from  the  screen  (x  >  0)  or  away  from  the  center  of  the  slab  (x  <  0). 
Examples  of  the  transverse  variation  of  leaky  waves  are  given  in  Figures 
3.3  and  3.6. 

The  launching  coefficient  contains  the  influence  of  the  excitation 
on  a  particular  mode.  In  this  respect  it  is  related  to  diffraction 
coefficients  computed  in  the  geometrical  theory  of  diffraction.  Vie  have 
computed  the  launching  coefficients  of  contributing  modes  for  several 

different  combinations  of  slab  parameters.  From  these  results  we  can 
gain  some  insight  into  the  factors  that  influence  the  launching  of 
guided  modes  in  the  slab  by  the  current  on  the  half  plane. 

g 

The  guided  wave  contributions  to  e ^  (x,y)  are  part  of  the  field 

diffracted  by  the  screen.  The  current  induced  on  the  half  plane  excites 

th>!  slab  modes  even  though  the  prir.ary  field  is  a  plane  wave.  For  the 

general  case  of  diffraction  by  a  half  plane,  diffraction  from  the  edge 

of  the  half  plane  has  the  most  influence  on  the  diffracted  field.  Thus, 

we  might  expect  that  the  edge  diffraction  effect  plays  the  major  role  in 

SDP 

the  launching  of  modes  in  the  slab.  Since  e^  (x,y)  is  heavily 
influenced  by  the  edge  diffraction  effect,  from  it  we  can  develop  an 
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zeros  of  an  analytic  function  in  a  specified  region  by  dividing  the 

region  into  a  number  of  subregions  and  then  performing  a  numerical 

integration  of  the  function  around  each  subregion.  Using  the  principle 

of  the  argument  (Hille,  1959b),  the  routine  determined  if  a  zero  was 

located  in  a  particular  subregion.  If  there  were  one  or  more  zeros  in  a 

subregion,  the  small  region  provided  a  good  approximate  starting  place 

for  applying  Muller's  method.  With  these  poles  located  it  was  a  simple 

+ 

matter  to  determine  if  they  were  swept  in  the  deformation  to  SDP  or 

_  p 

SDP  and  then  to  compute  any  resulting  contribution  to  e^  (x,y). 


3.3  Modal  Field  Components 

With  the  contributing  SWP's  and  LWP's  located,  we  can  compute  the 
launching  coefficient,  transverse  modal  function,  and  propagation  factor 
for  a  particular  pole. 

The  propagation  factor  (in  eq.  2.69)  for  a  surface  wave  shows 
propagation  in  the  positive  x  direction  for  x  >  0  and  in  the  negative  x 
direction  for  x  <  r  This  is  expected  since  we  expect  energy  to  be 

radiated  away  from  the  edge  of  the  half  plane.  If  the  slab  is  lossy  the 
propagation  factor  shows  some  decay  due  to  the  fact  that  moves  off 
the  vertical  line  through  ±u/2  so  that  sina^  becomes  complex.  For  a 
leaky  wave,  the  propagation  factor  shows  that  energy  is  being  radiated 
away  from  the  edge  of  the  half  plane,  but  it  also  incorporates  a 
significant  amount  of  decay  due  to  the  imaginary  part  of  sina^. 
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C(k  )  = 

X 

B2co6( B^t)  + 

j  8j  sin(  62  t)  =  0,  odd 

modes,  x>0. 

(3.5) 

C(k  )  = 

X 

B2cos[  &2  (-j)j 

+  j  Bj  sin[  &2  1  *  0, 

odd  modes ,  x<0. 

(3.6) 

C(k  )  = 

X 

8jC06[  e2  (-j)] 

+  j  828in[  BjOj)]  =  0, 

even  modes,  x<0. 

(3.7) 

By  using  top  sheet  values  of  B^  ,  we  can  find  the  SWP's,  and  by  using 
bottom  sheet  values  of  8^  ,  we  can  find  the  LWP's.  For  a  lossless  slab 
contributing  SWP's  lie  on  the  segment  (k^.k^)  for  x  >  0  or  on  the 
segment  (-k^-k^)  for  x  <  0  (Coblin  1983).  These  segments  map  into  the 
a  plane  as  segments  with  u  =  ±n/2  and  v  t  or*  v  e  (v_k2,0)  for 

x  >  0  or  x  <  0,  respectively,  where  v^  is  the  imaginary  part  of  the 
mapping  of  k^  into  the  a  plane  (Figures  2.3  and  2. A).  If  the  slab  is 
lossy,  the  SWP's  will  6till  lie  on  or  near  these  segments.  Contributing 
LWP's  lie  in  the  mapping  of  the  fourth  quadrant  of  the  bottom  sheet,  , 
for  x  >  0  and  in  the  mapping  of  the  second  quadrant  of  the  bottom  sheet, 
B^,  for  x  <0  as  shown  in  Figure  2.3.  According  to  Barone  (1956)  there 
are  an  infinite  number  of  LWP's,  but  only  a  finite  number  contribute. 
This  is  because  the  SDP  is  asymptotic  to  the  boundary  between  the 

+ 

mappings  of  the  top  and  bottom  sheets  (Figure  2.4).  Since  neither  SDP 
nor  SDP  sweeps  the  regions  B^  or  B2,  respectively,  at  infinity,  not  all 
the  LWP's  are  swept,  and  only  a  finite  number  contribute. 

The  SWP's  and  LWP's  were  found  using  a  numerical  search  routine 
developed  by  Singaraju,  Giri,  and  Baum  (1976).  This  routine  found  the 
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Under  the  screen  the  SWP's  and  LWP's  are  zeros  of 

C'(kx)  =  (61  -  62)e"j62t  -  (Bj  ♦  B2)e+j62t,  x  >  0.  (3.3) 


In  this  case  C'Ck^)  is,  within  a  factor,  simply  the  numerator  of 

G(k  ,0).  Thus,  the  poles  of  G/G  for  x  >  0  are  the  same  as  the  zeros  of 
x  + 

G(kx,0),  validating  our  claim  that  the  slab  shows  the  charac teristics  of 
a  grounded  dielectric  slab  waveguide  for  x  >  0.  For  x  <  0  the  SWP's  and 
LWP's  are  the  zeros  of 


C '  (k  )  =  (8  -  8  )2e  ^  B2 t  -  (  B.  +  B  )2e  +  ^B2t 

x  1  2  12 


s  ■^62(2)  -  (B,  +  B„)e+:’B2(2)  • 


=  (  Bj  -  B2)e  J  2V2  -  (  B1  +  B2! 


|(  8.  -  6,)e~-’62(2)  +  (8,  +  B,)e+:’B2(2)  ,  x  <  0. 

i,  1  2  l  l  ) 


(3.4) 


In  this  case  C'(k  )  is,  within  a  factor,  the  denominator  of  G(k  ,0). 

X  x 

Thus,  the  poles  of  G/G+  for  x  <  0  are  the  same  as  the  pole6  of  Gik^.O), 
that  is,  the  slab  shows  the  characteristics  of  a  dielectric  slab 
waveguide  in  space  for  this  region. 

From  the  expression  for  C'(k  )  we  can  develop  a  set  of 
characteristic  equations  for  each  case  (i.e.,  modes  odd  with  respect  to 
the  half  plane  for  x  >  0,  modes  odd  with  respect  to  the  center  of  the 
slab  for  x  <  0,  and  modes  even  with  respect  to  the  center  of  the  slab 
for  x  <  0).  These  characteristic  equations  can  be  written  as 


igure  3.1.  Generic  field  distributions  for  dielectric  slab  waveguide 
(a)  odd  modes  for  open  slab,  (b)  even  modes  for  open  slab 
(c)  "grounded-slab"  mode  under  p.e.c. 
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- EQUiPHASE  PLANES  - EQUIPHASE  PLANES 

- EQUI  AMPLITUDE  PLANES - EOUIAMPLITUDE  PLANES 


Figure  3.2.  Surface  and  leaky  waves:  (a)  surface  wave,  (b)  leaky  wave, 
(c)  region  of  contribution  for  a  leaky  wave. 
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Figure  3.8.  Representat i ve  launching  coefficients  for  a  slab 

with  i  1  =  10* q,  t*  =  Xq/ 2  showing  nonzero  behavior 
»  0°  and  "  =  180°. 


at  n 


BIBLIOGRAPHY 


Barone,  S.  (1956),  "Leaky  Wave  Contributions  to  the  Field  of  a  Line 

Source  above  a  Dielectric  Slab",  Microwave  Res.  Inst..,  Polytechnic 
Inst,  of  Brooklyn,  Report  R-532-56,  PIB-462,  November  1956. 

Coblin,  R.D.  (1983),  "Scattering  of  an  Electromagnetic  Plane  Wave  from  a 
Perfectly  Electrically  Conducting  Half-plane  in  the  Proximity  of 
Planar  Media  Discontinuities",  Ph.  D.  dissertation.  University  of 
Mississippi,  February  1983. 

Felsen,  L.B. ;  N.  Marcuvitz  (1973a),  Radiation  and  Scattering  of  Waves. 
Englewood  Cliffs,  N.J.:  Prentice-Hall,  pp.  278-284. 

Felsen,  L.B.;  N.  Marcuvitz  (1973b),  op.  cit.,  pp.  377-399. 

Felsen,  L.B.;  N.  Marcuvitz  (1973c),  op.  cit.,  pp.  399-410. 

Harrington,  R.F.  (1961),  Time-Harmonic  Electromagnetic  Fields.  New  York: 
McGraw-Hill. 

Heitman,  W.G. ;  P.M.  van  den  Berg  (1975),  "Diffraction  of  Electromagnetic 
Waves  by  a  Semi-infinite  Screen  in  Layered  Medium",  Can.  J.  Phys . , 
53(14),  pp.  1305-1317. 

Hessel,  A.  (1969a),  "General  Characteristics  of  Traveling-wave 

Antennas",  in  Collin,  R.E. ;  F.J.  Zucker ,  Antenna  Theory.  New  York: 
McGraw-Hill,  1969,  pp.  151-258. 

Hessel,  A.  (1969b),  op.  cit.,  pp  166-181. 

Hessel,  A.  (1969c),  op.  cit.,  pp  165-180. 

Hille,  E.  (1959a),  Analytic  Function  Theory,  Vol.  1.  New  York:  Chelsea, 
pp.  203-204. 

Hille,  E.  (1959b),  op.  cit.,  p.  253. 

Mittra,  R.  ;  S.W.  Lee  (1971a),  Analytic  Techniques  in  the  Theory  of 
Guided  Waves.  New  York:  Macmillan,  pp.  73-153. 

Mittra,  R.  and  S.W.  Lee  (1971b),  op.  cit.,  pp.  23-29. 


59 


60 


Singaraju  B.D.;  D.V.  Giri;  C.E.  Baum  (1976),  "Further  Developments  in 
the  Application  of  Contour  Integration  to  the  Evaluation  of  the 
Zeros  of  Analytic  Functions  and  Relevant  Computer  Programs",  Air 
Force  Weapons  Lab.,  Mathematical  Notes,  Note  42,  March  1976. 

Tamir ,  T. ;  L.B.  Felsen  (1965a),  "On  Lateral  Waves  in  Slab  Configurations 
and  Their  Relation  to  Other  Wave  Types",  IEEE  Trans.  Antennas 
Propagat.,  AP-13,  pp.  410-422. 

Tamir,  T. ;  L.B.  Felsen  (1965b),  op.  cit.,  pp.  419-422. 


DATE 

FILMED 


t 


